BI-TP 2013/10 



o 

Oh' 
< 

O 
m 



Oh- 
I . 

CD 



> 

O 
O 
00 

^' 

O 
m 



X 
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We evaluate the free energy of the Yang- Mills theory using the Gribov quantization that copes with 
non-perturbative resunimation. The magnetic scale is automatically incorporated in the framework 
and we find it efficient to stabilize the perturbative expansion of the free energy. In the range of 
the temperature T = Tc^2Tc major uncertainty in our results comes from the non-perturbative 
running coupling that is adopted from the lattice simulation, while the convergence above 2Tc is 
impressively robust. 
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Introduction — Stimulated by the exciting develop- 
ments from the Relativistic Heavy Ion Collider (RHIC) 
and the Large Hadron Collider (LHC) experiments, ther- 
modynamics of the quark-gluon plasma, especially the 
pressure p or the free energy density / ~ — p, is of cru- 
cial interest. Due to asymptotic freedom [1], perturba- 
tion theory was expected to work for QCD thermody- 
namics at high temperatures. Weak-coupling expansion 
of the QCD free energy has been accomplished up to 
g^logg, which is the highest order in the perturbative 
approach, with g being the running coupling (see [2] for 
reviews). Unexpectedly, the resulting weak-coupling se- 
ries shows poor convergence in the intermediate tempera- 
ture regime, i.c.T = 2T^^ATc with T^ ~ 170 Me V being 
the pseudo-critical temperature for the QCD phase tran- 
sition, which is probed in the RHIC and the LHC experi- 
ments. Resummation must be carried out to incorporate 
the contributions from the electric scale gT [3-6] . 

The non-convergence of QCD thermodynamics is at- 
tributed to the infrared (IR) sector of QCD which is gov- 
erned by soft-scale gluons [2]. Therefore the pure Yang- 
Mills (YM) theory provides the simplest test bench to ac- 
cess the IR problem. While the resummed perturbation 
theory in the electric sector resulted in improved conver- 
gence of the YM free energy down to T ^ 3 Tc, significant 
deviations from the lattice data were observed as further 
lowering T toward T^, which could be taken as a signal 
of the onset of non-perturbative effects. 

It was discovered in the early 1980s that due to 
the absence of an IR cutoff by screening in the mag- 
netic sector, perturbative expansion of finite-T YM the- 
ory breaks down at a fundamental level at the magnetic 
scale g^T. This IR catastrophe is the so-called Linde 
problem [7]. As a result, perturbation theory works up 
to a certain order only depending on the quantity in con- 
sideration. For the free energy IR divergences from the 
magnetic scale start entering at four loops correspond- 
ing to g^ order, which makes the highest accessible order 
g^Xogg as aforementioned. In the RHIC and LHC tem- 
perature regime the running coupling is g ~ £*(!), there- 
fore the missing contributions from the non-perturbative 



magnetic scale should be as significant as the ones from 
the perturbative electric scale. This fact provides us with 
an account for the aforementioned deviation in the free 
energy of resummed perturbation theory from the lat- 
tice data. There are various attempts to incorporate the 
magnetic contributions, but the non-perturbative nature 
is inherent in confining properties of dimensionally re- 
duced YM theory at high temperature [8] . We must thus 
augment the resummed perturbation theory with con- 
finement mechanism even when we deal with a deconfined 
state of matter. 

The key ingredients in constructing the YM thermo- 
dynamics are the correlation functions of the gluons and 
the ghosts. Gauge fixing is conveniently done through 
the Faddeev- Popov (FP) procedure [9], and the bare 
gluon and ghost propagators show an IR pole propor- 
tional to 1/p^ with p = {potP) being the Minkowskian 
four-momentum. The corresponding on-shell dispersion 
relation is po = |p|, which is modified in the full prop- 
agators including interaction effects; e.g. a plasmon ap- 
pears in the hard-thermal-loop resummation at the elec- 
tric scale. 

In order to have a better understanding in the 
intermediate temperature regime, the YM correlation 
functions must be improved non-perturbatively includ- 
ing the magnetic sector. The last decade has witnessed 
tremendous developments of high-precision lattice mea- 
surements of the YM correlation functions at both zero 
and finite temperatures [10]. There has been a consider- 
able progress in the study of Yang-Mills correlation func- 
tions also by functional methods in e.g. [13], and the 
results are in good agreement with the lattice simulation 
(see [14] and references therein). The obtained gluon and 
ghost propagators show an IR suppression and enhance- 
ment, respectively, as compared to the FP case. These 
results heuristically encompass desirable features of con- 
finement: The IR-suppressed gluon propagator indicates 
gluon confinement at large distance, and the IR-enhanced 
ghost is responsible for confinement. It has also been 
demonstrated that the balance in the IR sector is re- 
sponsible for the deconfinement phase transition [11, 12]. 



Gribov pointed out more than three decades ago 
that there still remains a residual gauge ambiguity in the 
IR sector, i.e. the so-called Gribov copy problem [15]. 
If the functional integration is dominated by the con- 
tributions near the horizon of the fundamental modular 
region, the ghost propagator is naturally enhanced in the 
IR momenta, and the resulting gluons are suppressed in 
turn. This is in line with the results from the lattice 
and functional methods. The corresponding on-shell dis- 
persion relation reads po = \/p^ + ttiq/p^ with ttig be- 
ing the Gribov mass which is solved by the condition of 
the horizon dominance. It was shown by Zwanziger in 
a phcnomenological way that a free gas of Gribov quasi- 
particles qualitatively captures the non-perturbative fea- 
tures of the lattice equation of state [16]. Subsequently, 
it was found that the Gribov mass is correctly propor- 
tional to the magnetic scale, i.e. tog ^ g^T in the limit 
of T -^ oo [17]. This finding is promising enough and 
confining mechanism can indeed resolve the Linde prob- 
lem. Many of the works on the Gribov copy problem are 
dedicated to the vacuum (see [18] for a review). System- 
atic attempts to extend the theory of confinement to the 
finite-T problems should deserve further investigations 
along the same line as a first try [17]. In this work we will 
pursue this possibility to make the perturbative evalua- 
tion of YM free energy stabilized by the non-perturbative 
mass emerging from the Gribov gauge-fixing procedure in 
which confinement is implemented, which we simply call 
the Gribov quantization. 

Formalism — In the Gribov quantization [15], the YM 
partition function in Euclidean space reads, 

Z= f VA{x)V{n)S{d-A) dct[-d-D{A)]e-^™ , (1) 
Jn 

in which the functional integration should be carried out 
within the Gribov region fi. 



n = { A: d-A^O , -d-D{A) > } , 



(2) 



whereby the Landau gauge is chosen. The restriction of 
the integration to the Gribov region is realized by insert- 
ing a function V{^) to the partition function (1), where 
V{Vi) represents the no-pole condition, so that the inte- 
gration never exceeds the Gribov horizon, i.e. 



v{n) = 6i[i-o-(o)] = 
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where 1 — a{P) is the inverse of the ghost dressing func- 
tion Zg{P) [see Eq. (11)]. The integration variable /3 is 
identified as the Gribov mass toq after a redefinition (see 
[18] for technical details of the Gribov quantization). 

Gap equation : The Gribov mass ttt-g is a Lagrange mul- 
tiplier that is determined by the variational principle, 
leading to the following gap equation as 
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where Nc represents the number of colors and d is the 
space-time dimensions; d = 4 — 2e. We note that P 
denotes the Euclidean four-momentum, P = ip,Pi) and 
P4 = 2mrT. The sum-integral is defined by 
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for convenience in our calculations using the dimensional 
regularization and the MS renormalization scheme. 

After carrying out the sum-integral and subtracting 
the UV divergence (see [19] for details) the gap equation 
becomes 
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where w± = ^Jp^ ± im^ and n^{x) = (e^'^ — 1)^"'^ is the 
Bose-Einstein distribution function. Analytical solutions 
are available in the limiting cases as. 
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in which it is evident that the magnetic scale emerges at 
high temperature. In what follows we will solve Eq. (6) 
numerically to derive tog as a function of T for a given 
renormalization scale fi (that will be fixed later). 

Gluon propagator : Gribov's giuon propagator in the 
Landau gauge reads. 
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which is regular as P — > featuring the IR suppression 
of confined gluons [15]. To proceed to any calculation 
involving Da{P), we need to specify the running coupling 
as a function of T. Since we are interested in the IR 
regime, the perturbative running is not quite appropriate 
and we have to use the full non-perturbative running. 
Instead of solving the non-perturbative resummation, in 
this work, we adopt the one measured on the lattice at 
finite T [20]. In [20] the running coupling is extracted 
from the large distance (IR) and the short distance (UV) 
behavior of the heavy-quark free energy leading to aiR 
and auv, respectively. 

Interestingly, the lattice measured coupling can be 
nicely fitted by a one-parameter form as 
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with c = 1.43 for the IR case and c = 2.97 for the UV case 
as shown together with the lattice data in Fig. 1. In the 
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FIG. 1. Lattice data of IR (upper) and UV (lower) running 
couplings at finite T from [20] and the fit results with a choice 
of cin Eq. (9). 
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FIG. 2. Ghost dressing functions at T = 200 MeV (by sohd 
curve) and T — 400 MeV (by dashed curve) as compared to 
the lattice data from [22]. 



intermediate temperature regime, the running coupling 
has a substantial scheme dependence and we therefore 
use the IR and the UV values to estimate theoretical 
uncertainties. 

Ghost propagator : The no-pole condition (not to go 
across the Gribov horizon) requires that the IR limit of 
the ghost propagator should be enhanced indicating the 
proximity to the Gribov horizon [15]. The ghost propa- 
gator in the Landau gauge thus reads, 
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where the ghost dressing function Zq = [1 — a{P)]^^ 
and cr(P) with gluon ladder diagrams taken into account 
turns out to be 



a{P) = N,g- 
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(11) 
Using the gap equation (4) one can confirm 1 — cr(0) = 

featuring the IR enhancement of the ghosts. We note 
that recent lattice data favor the so-called decoupling so- 
lution for which the ghost dressing function stays finite 
at P — >■ 0, which can be incorporated in a refined Gribov- 
Zwanziger approach [21]. In this study, however, we keep 
using the original form as above, so that we do not have 
to introduce any more adjustable parameters. Different 
behaviors in the deep IR region hardly affect bulk ther- 
modynamics. 

At T = this integration can be performed with the 
dimensional regularization, which yields, 
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where mco = '71g(T' = 0), and we find that a choice 
of /i = 1.69 GeV (in view of Eq. (7) with g = 3.13 ad- 



justed by hand) can reproduce the T = lattice data of 
Zq [22, 23]. Amazingly, with the same fj, fixed at T = 0, 
our finite-T numerical results from Eq. (11) [with g{T) 
from Eq. (9)] are rather insensitive to T, which is per- 
fectly consistent with the lattice simulation. We note 
that in the direct evaluation of Eq. (11) we imposed a 
three-momentum cutoff A = 1.25/^ so that the T = nu- 
merical results are matched with Eq. (12). Fig. 2 shows 
the numerical results at T = 200 MeV and T = 400 MeV 
when g{T) runs with c = 2.97. If we use c = 1.43, for the 
T = 200 MeV case, cT is too close to T^. As long as T 
is greater than Tc, however, the behavior of Zq is robust 
and not contaminated by the uncertainty in g(T). 

Now that we have confirmed that the resummed 
ghost dressing function in our approach is almost T- 
independent as observed in the lattice simulation, we can 
safely use Eq. (12) also for the finite-T calculation to re- 
duce the computational cost. 

Quasi-particle Approximation — With the gluon 
and ghost propagators prepared, we can calculate the 
free energy. In the 2-particle irreducible (2PI) formalism 
the effective action can be expressed in terms of the full 
propagator G and full self-energy 11 as 
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where r2[G] represents the sum of all 2PI diagrams. 
In a general quasi-particle approximation we keep only 
the first term, which should be a reasonable estimate; 
the variational principle to derive G can be formulated 
to minimize the sensitivity to higher-order corrections, 
which is often called the optimized perturbation theory 
(see [24]). 

Gluon contribution : Using the gluon propagator (8) the 
gluon contribution to the free energy reads. 
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FIG. 3. Gribov mass as a function of T. The band shows the 
uncertainty arising from different schemes of lattice-measured 
aa{T) (see Fig. 1). The upper (solid) and lower (dashed) 
bounds correspond to the IR and UV couplings, respectively. 



which can be evaluated straightforwardly in the dimen- 
sional regularization and the results are expressed as 



^trlnD^' = {N^-l)^ 
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Ghost contribution : The ghost contribution to the free 
energy reads, 



- tvlnD;' = -{Nl - 1)^^ ln{p2 [ i _ ^(p) ] } 
^ ln[l-a(P)] 



= {N^ - 1) 



^2t4 



(16) 



with Eq. (12) substituted for l — a{P) as justified numer- 
ically. It is technically daunting to do the sum-integral of 
Eq. (16) analytically, and we will resort to the numerical 
calculation, in which the UV divergence is subtracted at 
a certain T sufficiently below Tc. 

Results and Discussions — With the renormalization 
scale ^ fixed as ^ — 1.69 GeV, we can numerically solve 
the gap equation (6) to obtain toq as a function of T as 
shown in Fig. 3. 

Solving the gluon contribution (15) and the ghost 
contribution (16) numerically with ttiq as determined in 
Fig. 3, we obtain the free energy (or pressure) of Gribov 
quasi-partieles scaled by the massless Stefan-Boltzmann 
limit -(iV2 _ 1)7^2^4/45 f^^ iV^ = 3 as shown in Fig. 4. 
The SU(3) Yang-Mills lattice data (blue dots) are taken 
from [25]. The red band shows the uncertainty arising 
from the lattice running couplings with the lower (solid) 
and the upper (dashed) bounds corresponding to the IR 
and the UV couplings, respectively, as we have seen in 
Fig. 9. 

It is worth noting that our calculations are done 
in a way consistent with the lattice- measured as{T), 
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FIG. 4. Free energy (pressure): the red band represents our 
numerical results with the solid and the dashed curves corre- 
sponding to those in Figs. 1 and 3. The blue dots represent 
lattice data from [25] , and the gray band represents three-loop 
HTLpt results from [6]. 



while the perturbative running coupling with an assumed 
renormalization scale used in conventional resummed 
perturbation theory is not compatible with the lattice 
as{T) even at very high T (see discussions in [20] for 
more details). Due to the ambiguity in defining as{T) 
on the lattice, we see a rather wide band at T < 2Tc. 
However, the uncertainty in the resulting free energy gets 
suppressed significantly at T > 2.5 Tc, which is highly 
non-trivial especially considering the fact that there is 
still a big variation between the IR and the UV lattice 
couplings at such temperatures. 

Comparing to the most recent estimate from three- 
loop hard-thermal-loop perturbation theory (HTLpt) [6] 
as shown by the gray band in Fig. 4, the uncertainty 
in our free energy is about 35% of that of the HTLpt 
at 2.5 Tc, and about 15% at 5Tc. We would like to 
point out, furthermore, that our free energy is consis- 
tent with the lattice data in the whole displayed tem- 
perature range. In particular the lattice data lie in our 
band even below 3.5 T^, where HTLpt suffers from poor 
convergence, and thereby, sharp rising of the free energy 
after the dcconfinemcnt phase transition which features 
a non-perturbative release of new degrees of freedom is 
realized somehow even without inclusion of the Polyakov 
loop. The free energies from the resummation schemes of 
[3, 4] have relatively smaller uncertainties than HTLpt, 
but they are not successful enough to recover the shape 
rising behavior below 2Tc^ where we consider that the 
resummation in the magnetic sector is crucial. 

In closing, we would like to stress that the resumma- 
tion in the magnetic sector is by itself set in through the 
Gribov quantization, though the quasi-particle approxi- 
mation may look like a one- loop integration. In fact, the 
construction of cr(P) explicitly involves the higher-order 
diagrams, and solving the gap equation is nothing but the 
variational calculation to carry out the non-perturbative 



resummation effectively. In this sense, it is neither an 
accident nor a magic that our results are such consistent 
with the lattice data, but rather a natural consequence 
from proper reformation of the magnetic contributions. 

Conclusions and Outlook — In this Letter we have 
explored a novel and systematic evaluation of Yang-Mills 
free energy using the Gribov quantization. The results, 
thanks to the improvement at the magnetic scale, show 
robust and stable behavior consistent with the lattice 
data, which evidently manifests profound importance of 
the non-perturbative gauge fixing even in the perturba- 
tive evaluation of thermodynamics at high temperature. 
There are intriguing future extensions: First, the 
explicit calculation of the magnetic screening mass and 
the spatial string tension would provide us with a deeper 
insight into the Lindc problem. Second, the Polyakov 
loop coupling would enable us to investigate the nature of 
the first-order phase transition of deconfinement. Third, 
the running coupling from the functional renormaliza- 
tion group approach would help with reducing the uncer- 
tainty band [26, 27]. Fourth, it would be desirable to sys- 
tematize the perturbation theory using the local Gribov- 
Zwanziger action [28]. Last but not least, the application 
of the Gribov quantization would have a profound impact 
on more general QCD problems, e.g. quarkyonic chiral 
spirals with and without magnetic fields [29], transverse 
dynamics of real-time evolution [30], and so on. 
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